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IN THE SPECIFICATION 



Please amend the specification as follows: 
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Please replace page 4 with the following: 

fracture simulators. Methods employing PL3D aim at -accurately 
take — taking into account geologic layers . One such program 
prcocntod ao PL3D , known commercially as GOHFER (GOHFER is 
believed to be a trademark of Stim-Lab and the Marathon Oil 
Company) , provides grid oriented hydraulic fracture replication 
capabilities. This grid oriented program, and its mode of 
operation, is seen in Figure 7. As the front of the fracture 
moves forward, calculations are made in which each individual 
grid is either "on" or "off" depending upon whether or not more 
than half of the individual grid is "covered" by the advancing 
fracture as it moves outward from the wellbore. If more than 
one-half of the grid element is covered, then the element is 
estimated to be fully active. The disadvantage of this system of 
estimating fracture growth is that it produces too much numerical 
noise at the fracture tip, and hence in the output data. 
Other PL3D methods of simulating fractures include the TerraFrac 
three dimensional fracturing simulator (TerraFrac is a trademark 
of the TerraTek Company) . This simulator operates as seen in 
Figure 8, using estimates that are based upon a method of a 
moving mesh. This method shows less noise than the GOHFER 
method, because it uses triangle 
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Please replace page 22 with the following: 

fracture width and fracture pressure on each active element. A complete 
description of the process of the propagation of a hydraulic fracture is thus 
obtained. 

Solutions of the multi-layer equilibrium equations are provided. A three- 
dimensional body is assumed. The theory also applies to the two-dimensional, cases 
(plane strain, plane stress, antiplane strain). The method provides an efficient 
way of determining the solution to the equilibrium equations: 

(Tgj+b^O (1) 

for a transversely isotropic elastic medium with a stress strain 
relationship given by: 

<?U =C W € ki (2) 
In Equation (1) and (2) the index notation is used. In this standard 

notation, a repeated literal index in any term of an expression implies 

summation. Therefore, in equation (1) 0 , -3 ij J means 

-CFjj j= <7 in +Cj22 +c7 /33 since il is assumed a 3D body. 

The comma preceding an index denotes partial differentiation with respect 

to that variable represented by that index.;- u i tj thus denotes -? d u i / -?r)x^ . 

The notation used in Equation (1) is consequently a shortcut for describing 
the following equations: 

d&xx ^ a xy d<T xz . 

dx dy dz 
do vx 3cr w dcr vz 
dx dy dz y 

dx dy dz 
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Please replace page 25 with the following: 
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case of isotropic layers, the system (5) has the general . solution : 



(7) 



Here d Jr and f Jr are constants that depend on the material constants 



of the layer, the CJT are the roots of the characteristic equation for the 
system of ordinary differential equations, and the Aj(fc)are free parameters 
of the solution that are determined by the forcing terms ~bir hi in (1) and 

"the interface conditions prescribed at the boundary between each of the 

layers (e.g. bonded, f rictionless , etc.). 

Substituting these displacement components into the stress strain law 
(2), we can obtain the corresponding stress components: & xx , G ^ , & u , &^ , 

& , and & , which can be expressed in the form: 




(8) 
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Please replace equation (9) with the following: 
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